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Abstract 

OO ' In a recently developed approximation technique for quantum field theory 

On ' 

\ the standard one-loop result is used as a seed for a recursive formula that gives 

a sequence of improved Gaussian approximations for the generating functional. In 
O \ this paper we work with the generic (f)^ + (j)^ model in d = dimensions. We compare 

the first, and simplest, approximation in the above sequence with the one-loop and 
two- loop approximations, as well as the exact results (calculated numericaly). 



o 

in 



> 



The central object in quantum field theory is the generating functional Z[J]. Func- 



CN \ tional derivatives of Z[J] with respect to the external fields J(x) give the Green's functions 

of the theory. The generating functional is determined from the (Euclidian) action S[(j)\ 
through the path integral 

00 

g ; Z[J] = I [d<P] e-(5M-/'^-^W<^W) . (1) 

Qh! The integration measure is, formaly, simply = HaiGiRd d(l>{x), where d is the dimension 

^ . of space-time. In this paper we will work with models in = dimensions. In d = 

functionals become functions, and the path integral reverts to a single definite integral 
^ ■ over the whole real line 

; Z{J) = [ d0e-(^(*)--^<^) . (2) 



Two further important objects are W{J) — the generator of connected diagrams (or free 
energy) 

Z( J) = Z(0) e-^(-^) , (3) 

and the quantum average of the field ip = {(p) = — jj W{J). 

In the Gaussian approximation, we Taylor expand the action in the path integral 
around some reference point (/)ref, and keep terms that are at most quadratic in — ^^ef- 
The integral in Q is now a Gaussian and we find 

Wg.UJ, M = S{M - J0ref + \ ln5"(0ref) " \ ^^^^^jf] "/^^ • (4) 

2 2 S"(0i.cf) 
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For this approximation to make sense, the integral must get its dominant contribution 
from the vicinity of the reference point (/ij-ef • The standard Gaussian approximation (loop 
expansion) corresponds to the choice 0ref = 0ciass (J), where 0ciass is the solution of the 
classical equation of motion 

S' = J. The classical solution is the maximum of the integrand in 
In a previous paper [Tj we expanded the integrand around the average field (p. As 
we have shown, although the classical solution gives the maximum of the integrand, 
expansion around ip gives a better approximation for the area under the curve. The 
Gaussian approximation around the average field ip is simply 



(5) 



To be able to calculate this in closed form we need to know (p{J), which is tantamount to 
knowing how to do the theory exactly, since Lp and its derivatives give all the connected 
Green's functions. The use of equation (jSJ comes about when one solves it iteratively. 
Using the definition of (p in terms of W, as well as equation Q we obtain the following 
iterative process 

d 

(6) 
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For the seed of this iteration we chose the classical field, i.e. </?o = 0ciass- In this way 
one obtains a sequence of points ipQ,ipi,ip>2i ■ ■ ■ or equivallently of approximations to the 
connected generating functional Wi, W2, W3, . . . given by Wn+i{J) = WQa.uss{J, fni^J))- In 
P] we have shown that this sequence gives better and better approximations and converges 
(though slowly) to the best Gaussian approximation ipoo- This is shown in Figure ^ Note 
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Figure 1: Plots of — y^o (dotted line), ip — ipi (dashed line), 'p — 'p)2 (thin line) and 'p — 'Pc 
(thick line) as functions of J. The action is given in (jH)) with couplings (73 = 0, (74 = 1. 



that Lpao is still not equal to the exact result Lp. The reason for this is obvious: We used the 
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Gaussian approximation VFcausa in defining our recursive relation, and there is no reason 
to expect that this converges to the exact resuh. 

So far we have seen that we can improve on the usual loop expansion. However, the 
sequence of improved Gaussian approximations converges very slowly. For this reason 
it is interesting to look at the first approximations in this sequence and compare them 
to standard approximation schemes. As we have seen, ipo = (;/)ciass, so Wi{J) is just the 
one-loop result. Our first new approximation is therefore W2{J) = VFGauss(<^i V'l)) where 
ipi = —-^Wi. The nice thing about this approximation is that it is only a bit more 
complicated than the one- loop result. From now on we will study this approximation, 
and compare it to one-loop and two-loop results, as well as to the exact results that 
have been calculated numericaly. In what follows we will designate W2 as the improved 
Gaussian approximation. 

For our comparison we have looked at the (full) n-particle Green's functions 
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-5(0) 



J d(j)e-'^('t>) 



(7) 



The model we considered was 



5(0) = ^0' + ^^730' + ^^740^ (8) 

All our calculations have been done for > ^g^ where the above action has a unique 
minimum at = 0. For the two-point function the one-loop approximation is given 
in Figure El In this plot we see the expected perturbative region in the vicinity of the 
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Figure 2: Absolute values of relative errors of the one- loop approximation to G2. Only 
the contours corresponding to |5G'2~^°°''| < 0.1 are shown. 
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origin. What is not immediately obvious is the meaning of the two regions sprouting off 
to large values of the coupling constants. Still, it is rather easy to give a simple hand- 
waving argument: The one-loop result for Gn is an (n — 1)^* order polynomial in the 
couplings, while the exact result is a relatively slowly varying, monotonous function (on 
the range of interest). As a consequence, SG]^^""^ can vanish on (at most) n curves in the 
5^3, (/4 plane. The sprouting regions in Figure |21 flank these two curves. A similar plot of 
5G\^^°°^ has one sprouting region. What is important to note is that the sprouting regions 
corresponding to different Green's function have nothing to do with each other. In fact, 
they are just a manifestation of the old saw that a stopped clock gives the correct time 
of day twice during each day. The same plot for the two-loop and improved Gaussian 
approximations is shown in Figure El As before, we again have central (perturbative) 
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Figure 3: Contour plot of < 0.1 (left), and I^G™™*^] < 0.1 (right) 

regions and a certain number of sprouting regions. These plots clearly show that the 
improved Gaussian approximation outperforms both the one-loop and two-loop results. 
Figure El gives the same comparison of one- loop, two- loop and improved results for the 
one-point Green's function. 

We have looked at all the Green's function from Gi to G^. In all cases the improved 
Gaussian approximation gives the best result, however the advantage becomes less marked 
when one looks at higher Green's function. 

As a conclusion, we have shown in this paper that even the simplest improved Gaussian 
approximation gives better agreement with exact results than both one-loop and two-loop 
approximations. At the same time, the computational cost of the improved Gaussian 
approximation is negligibly greater than that of the one- loop result, and significantly 
smaller than of the two-loop result. In our previous paper we looked at the benefits of the 
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Figure 4: I^G^'"""! < 0.1 (top left), |5G'^'°°P| < 0.1 (top right), |5G'™™''| < 0.1 
(bottom) 

improved Gaussian approximation both from the analytic and numerical (Monte Carlo) 
sides. We are currently working |2I on extending those results to interesting models in 
d>l. 
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